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ABSTRACT 

Many of the current anomalies reported in the Wilkinson Microwave Anisotropy Probe 
(WMAP) l- vear data disappe ar after 'correcting' for the best-fit embedded Bianchi type Vllh 



component ( [laff^e et alj[2QQ5D , albeit assuming no dark energy component. We investigate the 



eff'ect of this Bianchi correction on the detections of non-G aussianity in the WMA P data that 
we previously made using directional spherical wavelets ( McEwen et al. 2005 a| ). We con- 
firm that the deviations from Gaussianity in the kurtosis of spherical Mexican hat wavelet 
coeflftcients are elir ninated once the data is corrected for the Bianchi component, as previ- 
ously discovered by Taff^e et al. ( 2005 ). This is due to the r eduction of the cold spot at Galactic 
coordinates (/, b) = (209°, -57°), which [Cruz et al. (2005) claim to be the sole source of non- 
Gaussianity introduced in the kurtosis. Our previous detections of non-Gaussianity observed 
in the skewness of spherical wavelet coefficients are not reduced by the Bianchi correction. 
Indeed, the most significant detection of non-Gaussianity made with the spherical real Morlet 
wavelet at a significant level of 98.4% remains (using a very conservative method to estimate 
the significance). Furthermore, we perform preliminary tests to determine if foregrounds or 
systematics are the source of this non-Gaussian signal, concluding that it is unlikely that these 
factors are responsible. We make our code to simulate Bianchi induced temperature fluctua- 
tions publicly available. 

Key words: cosmic microwave background - methods: data analysis - methods: numerical 
- cosmology: model 



1 INTRODUCTION 

Recent measurements of the cosmic microwave background 
(CMB) anisotropics made by the Wilkinson Microwave Anisotropy 
Probe (WMAP) provide full- sky data of unprecedented precision 
on which to test the standard cosmological model. One of the most 
important and topical assumptions of the standard model currently 
under examination is that of the statistics of the primordial fluc- 
tuations that give rise to the anisotropics of the CMB. Recently, 
the assumption of Gaussianity has been questioned with many 
works highlighting deviations from Gaussianity in the WMAP 1- 
year data. 

A wide range of Gaussianity analyses have been performed on 
the WMAP 1-year data, calculati ng measures such as the bispec- 
trum and Minkowski functionals (iKomatsu et al||2003|; M agueijo 
& Me deiros |2004|; [Land & Magueijol |2005a|), the genus ( Cofley 
& Gott |2003[ ^riksen et al.| |2004), correlation functions (Gaz - 
tanaga & Wagg |2003| ; ^riksen et al] |2005|; [Tojeiro et al.| |2o"o6|) 
low-multipole alignment statistics ( |de Oliveira-Costa et alj [2004 ; 



Copi et al, |2004 12006|; Schwarz et alj|2004; 


|Slosar&Seljak||2004; 


Weeks |2004|; |Land & Magueijc| 




ill 


[Bielewicz et aljp005), 



etal. |2005[ ), local curvature ( [Hansen et al.| |2004); pabelTa etal^ 
2005), the higher cr iticism statistic ([Cayon et ar _ 2005[), hot and 



cold spot statistics ( [Larson & Wandelj Eool [2005[; ICruz et al 
2005) and w avelet coeflicient statistics "[Vielvaetal.[[2003f Mukher- 
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jee & Wang [2004 [McEwen et alj [2005a[ ). Some statistics show 
consistency with Gaussianity, whereas others provide some evi- 
dence for a non-Gaussian signal and/or an asymmetry between the 
northern and southern Galactic hemispheres. Although these de- 
tections may simply highlight unremoved foreground contamina- 
tion or other systematics in the WMAP data, it is important to also 
consider non-standard cosmological models that could give rise to 
non-Gaussianity. 

One such alternative is that the universe has a small universal 
shear and rotation - these are the so-called Bianchi models. Re- 
laxing the assumption of isotropy about each point yields more 
complicated solutions to Einstein's field equations that contain 
the Friedmann-Robertson- Walker metric as a special case. B arrow 
et al. ( [1985[ ) derive the induced CMB temperature fluctuations that 
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result in the Bianchi models, however they do not include any dark 
energy component as it was not considered plausible at the time. 
There is thus a need for new derivations of solutions to the Bianchi 
models in a more modern setting. Nevertheless, the indu ced CMB 
temperature fluctuations derived by Barrow et al. (1985) provide a 



good phenomenological setting in which to examine and raise the 
awareness of more exotic cosmological models. 

Bianchi type Vllh models have previously been compared 
both to the Cosmic Backgr ound Explorer-Diff erential Microwave 
Radi ometer (COBE-DMR) (iKogut et al. || l997[) and WMAP (J affe 



et al. ]2005| ; henceforth referred to as J05) data to place limits on the 
global rotation and shear of the universe. Moreover, JOS find a sta- 
tistically significant correlation between one of the Bianchi Vllh 
models and the WMAP internal linear combination (ILC) map. 
They then 'correct' the ILC map using the best-fit Bianchi tem- 
plate and, remarkably, find that many of the reported anom a lies in 
the WMAP data disappear. More recently Land & Magueijc] ( 2006| ) 
perform a modified template fitting technique and, although they 
do not report a statistically significant template fit, their corrected 
WMAP data is also free of large scale anomalies. 

In this paper we are interested to determine if our previ- 
ous detection s of non-Gaussianity m ade using directional spheri- 
cal wavelets (McEwen et al. 2005a) are also eliminated when the 



WMAP data is corrected for the best-fit Bianchi Vllh template de- 
termined by J05. In section ^ the best-fit Bianchi template embed- 
ded in the WMAP data is described and used to correct the data, 
before a brief review of the analysis procedure is given. Results are 
presented and discussed in section ^. Concluding remarks are made 
in section 0. 



2 NON-GAUSSIANITY ANALYSIS 

We recently made significant d etections of non-Gaus sianity using 
directional spherical wavelets (McEwen etal. 20052 ) and are in- 



terested to see if these detections disappear when the data are cor- 
rected for an embedded Bianchi component. The best-fit Bianchi 
template and the correction of the data is described in this section, 
before we review the a nalysis procedure. We essentially repeat the 
analysis performed by McEwen et al. (20052) for the Bianchi cor- 



rected maps, hence we do not describe the analysis procedure in 
any detail here but give only a very brief overview. 



2.1 Bianchi Vllh template 

We have implemented simulations to compute the Bianchi-induced 
temperature fluctuations, concentrating on the Bianchi type Vllh 
models, which include the types I, V and VIIp as special cases ( Bar- 
row et al. 



1985[ ).^ Note that the Bianchi type Vllh models apply to 



open or flat universes only. In appendix Q we describe the equations 
implemented in our simulation; in particular, we give the analytic 
forms for directly computing Bianchi-induced temperature fluctu- 
ations in both real and harmonic space in sufficient detail to repro- 
duce our simulated maps. The angular power spectrum of a typical 
Bianchi-induced temperature fluctuation map is illustrated in Fig. |l| 
(note that the Bianchi maps are deterministic and anisotropic, hence 
they are not fully described by their power spectrum). An example 



^ Our code to produce simulations of Bianchi type Vll-induced tem- 
perature fluctuations may be found at: h ttp : //www . mrao . cam . ac . uk/ 
~jdm57/ 




Figure 1. Angular power spectrum of the Bianchi-induced temperature fluc- 
tuations. The particular spectrum shown is for the best-fit Bianchi template 
matched to the WMAP data. Notice that the majority of the power is con- 
tained in multipoles below ^ ~ 20. 



of the swirl pattern typical of Bianchi-induced temperature fluctu- 
ations may be seen in Fig. | (a) (this is in fact the map that has 
the power spectrum shown in Fig. y). Notice that the Bianchi maps 
have a particularly low band-limit, both globally and azimuthally 
(i.e. in both i and m in spherical harmonic space; indeed, only those 
harmonic coefficients with m-±\ are non-zero). 

The best-fit Bianchi template that we use to correct the 
WMAP data is simulated with the parameters determined by J05 
using the latest shear and vorticity estimates (J05; private commu- 
nication). This map is illustrated in Fig, ^ (d). In our previous non- 
Gaussianity analysis ([McEwen et al.||2005a| ) we considered the co- 
added WMAP map (Komatsu et al. p003[ ). However, the template 



fitting technique performed by J05 is only straightforward when 
considering full-sky coverage. The Bianchi template is therefore 
matched to the full- sky ILC map since the co- added WMAP map 
requires a galactic cut. Nevertheless, we consider both the ILC 
and CO- added WMAP map hereafter, using the Bianchi template 
matched to the ILC map to correct both maps. The Bianchi tem- 
plate and the original and corrected WMAP maps that we consider 
are illustrated in Fig. ^ 



2.2 Procedure 

Wavelet analysis is an ideal tool for searching for non-Gaussianity 
since it allows one to resolve signal components in both scale 
and space. The wavelet transform is a linear operation, hence the 
wavelet coefficients of a Gaussian map will also follow a Gaussian 
distribution. One may therefore probe a signal for non-Gaussianity 
simply by looking for deviations from Gaussianity in the distribu- 
tion of the wavelet coefficients. 

To perform a wavelet analysis of full- sky CMB maps we apply 
our f ast continuous spherical wavelet transform (CSWT) (M cEwen 
et al. [2005b| ), which is based on the spherical wavelet transform 
developed by Antoine, Vandergheynst and colleagues (Antoine & 



2005J 



lerghe 



Vander gheynst |1998| , |1999| ; [Antoine et al.|[2002^ |2004| - ]Wiaux et al 



and the fast spherical convolution developed by Wandelt 
& Gorski ( [200 1[ ). In particular, we use the symmetric and elliptical 
Mexican hat and real Morlet spherical wavelets at the scales de- 
fined in Table |lj The elliptical Mexican hat and real Morlet spher- 
ical wavelets are directional and so allow one to probe oriented 
structure in the data. For the directional wavelets we consider five 
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(a) Best-fit Bianchi template (scaled by 
four) rotated to the Galactic centre for il- 
lustration 



(b) ILC map 



(c) WMAP co-added map (masked) 




(d) Best-fit Bianchi template (scaled by 
four) 



(e) Bianchi corrected ILC map 



(f) Bianchi corrected WMAP co-added 
map (masked) 



Figure 2. Bianchi template and CMB data maps (in mK). The Bianchi maps are scaled by a factor of four so that the structure may be observed. The KpO 
mask has been applied to the co-added WMAP maps. 



evenly spaced azimuthal orientations between [0,;7r). We look for 
deviations from zero in the skewness and excess kurtosis of spher- 
ical wavelet coefficients to signal the presence of non-Gaussianity. 
To provide confidence bounds on any detections made, 1000 Monte 
Carlo simulations are performed on Gaussian CMB realisations 
produced from the theoretical power spectrum fitted by the WMAP 
team.^ The ILC map, the foreground corrected WMAP maps re- 
quired to create the co-added map, the masks and power spectrum 
may all be downloaded from the Legacy Archive for Microwave 
Background Data Analysis (LAMBDA) website^. 



3 RESULTS AND DISCUSSION 

We examine the skewness and excess kurtosis of spherical wavelet 
coefficients of the original and Bianchi corrected WMAP data to 
search for deviations from Gaussianity. Raw statistics with corre- 
sponding confidence regions are presented and discussed first, be- 
fore we consider the statistical significance of detections of non- 
Gaussianity in more detail. Localised regions that are the most 
likely sources of non-Gaussianity are then examined. Finally, we 
investigate the possibility of foreground contamination and system- 
atic s. 



3.1 Wavelet coefficient statistics 

For a given wavelet, the skewness and kurtosis of wavelet coeffi- 
cients is calculated for each scale and orientation, for each of the 
data maps considered. These statistics are displayed in Fig. ^ with 
confidence intervals constructed from the Monte Carlo simulations 
also shown. For directional wavelets, only the orientations corre- 
sponding to the maximum deviations from Gaussianity are shown. 



served bv Vielva et al. 


McEwen et al. 


(2005a|) 



( |2005a| ) in the kurtosis of the Mexican hat wavelet 



coefficients is reduced when the data are corrected for the Bianchi 
template, confirming the results of J05. However, it appears that 
a new non-Gaussian signal may be detected in the kurtosis of the 
symmetric Mexican hat wavelet coefficients on scale - 450' 
and in the kurtosis of the elliptical Mexican hat wavelet coeffi- 
cients on scale - 600'. These new candidate detections are 
investigated further in the next section. Interestingly, the skewness 
detections that we previously made are not mitigated when making 
the Bianchi correction - the highly significant detection of non- 
Gaussianity previously made with the real Morlet wavelet remains. 

It is also interesting to note that the co-added WMAP and ILC 
maps both exhibit similar statistics, suggesting it is appropriate to 
use the Bianchi template fitted to the ILC map to correct the co- 
added WMAP map. 



^ We use the theoretical power spectrum of the Lambda Cold Dark Matter 
(ACDM) model which best fits the WMAP, Cosmic Background Imager 
(CBI) and Arcminute Cosmology Bolometer Array Receiver (ACBAR) 
CMB ^m, 



ittp : //cmbdata . gsf c . nasa . gov/ 



3.2 Statistical significance of detections 

We examine the statistical significance of deviations from Gaus- 
sianity in more detail. Our first approach is to examine the distribu- 
tion of the statistics that show the most significant deviation from 
Gaussianity in the uncorrected maps, in order to associate signif- 
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Table 1. Wavelet scales considered in the non-Gaussianity analysis. The overall size on the sky ^\ for a given scale are the same 
for both the Mexican hat and real Morlet wavelets. The size on the sky of the internal structure of the real Morlet wavelet ^2 is 
also quoted. 



Scale 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


Dilation a 




100' 


150' 


200' 


250' 


300' 


350' 


400' 


450' 


500' 


550' 


600' 


Size on sky ^\ 


14V 


282' 


424' 


565' 


706' 


847' 


988' 


1130' 


1270' 


1410' 


1550' 


1690' 


Size on sky ^2 


i5.r 


31.4' 


47.1' 


62.8' 


78.5' 


94.2' 


110' 


126' 


141' 


157' 


173' 


188' 






(a) Skewness - Mexican hat e = 0.00 



(b) Skewness 
7 = 72° 



Mexican hat e = 0.95; 



(c) Skewness 
7 = 72° 



real Morlet A: = (10,0)^ 




(d) Kurtosis - Mexican hat e = 0.00 



(e) Kurtosis - Mexican hat e = 0.95; (f) Kurtosis - real Morlet A: = (10,0)^; 

7 = 108° 7 = 72° 



Figure 3. Spherical wavelet coefficient statistics for each wavelet and map. Confidence regions obtained from 1000 Monte Carlo simulations are shown for 
68% (red/dark-grey), 95% (orange/grey) and 99% (yellow/light-grey) levels, as is the mean (solid white line). Statistics corresponding to the following maps 
are plotted: WMAP combined map (solid, blue, squares); ILC map (solid, green, circles); Bianchi corrected WMAP combined map (dashed, blue, triangles); 
Bianchi corrected ILC map (dashed, green, diamonds). Only the orientations corresponding to the most significant deviations from Gaussianity are shown for 
the Mexican hat e = 0.95 and real Morlet wavelet cases. 



icance levels with the detections. Our second approach is to per- 
form tests on the statistics computed with each type of spheri- 
cal wavelet. This approach considers all statistics in aggregate, and 
infers a significance level for deviations from Gaussianity in the 
entire set of test statistics. 

Histograms constructed from the Monte Carlo simulations for 
those test statistics corresponding to the most significant deviations 
from Gaussianity are shown in Fig. 0. The measured statistic of 
each map considered is also shown on the plots, with the number 
of standard deviations these observations deviate from the mean. 
Notice that the deviation from the mean of the kurtosis statistics is 
considerably reduced in the Bianchi corrected maps, whereas the 
deviation for the skewness statistics is not significantly aflfected. 

Next we construct significance measures for each of the most 
significant detections of non-Gaussianity. For each wavelet, we de- 



termine the probability that any single statistic (either skewness or 
kurtosis) in the Monte Carlo simulations deviates by an equiva- 
lent or greater amount than the test statistic under examination. If 
any skewness or kurtosis statistic^ calculated from the simulated 
Gaussian map - on any scale or orientation - deviates more than 
the maximum deviation observed in the data map for that wavelet, 
then the map is flagged as exhibiting a more significant deviation. 
This technique is an extremely conservative means of construct- 
ing significance levels for the observed test statistics. We use the 
number of standard deviations to characterise the deviation of the 
detections, rather that the exact probability given by the simula- 



^ Although we recognise the distinction between skewness and kurtosis, 
there is no reason to partition the set of test statistics into skewness and 
kurtosis subsets. The full set of test statistics must be considered. 
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V ^AP 
(-5.61a) 

ILC 
(-4.69a) 

VHAP-BIANCm 
(-5.66a) 


1 




ILC-BIANCHI 
(-4.74a) 

J 





(a) Skewness - Mexican hat e = 0.00; 
«2 = 100' 



(b) Skewness - Mexican hat e = 0.95; as ■ 
150^7 = 72° 



(c) Skewness - real Morlet k = (10,0)^; 
an = 550^7 = 72° 





(d) Kurtosis - Mexican hat e = 0.00, 
ae = 300' 



(e) Kurtosis - Mexican hat e 
aio = 500'; 7 = 108° 



0.95; (f) Kurtosis - real Morlet A: = (10,0)^; 

an = 550'; 7 = 72° 



Figure 4. Histograms of spherical wavelet coefficient statistic obtained from 1000 Monte Carlo simulations. The mean is shown by the thin dashed black 
vertical line. Observed statistics corresponding to the following maps are also plotted: WMAP combined map (solid, blue, square); ILC map (solid, green, 
circle); Bianchi corrected WMAP combined map (dashed, blue, triangle); Bianchi corrected ILC map (dashed, green, diamond). The number of standard 
deviations these observations deviate from the mean is also displayed on each plot. Only those scales and orientations corresponding to the most significant 
deviations from Gaussianity are shown for each wavelet. 



tions, since for many of the statistics we consider no simulations 
exhibit as great a deviation on the particular scale and orientation. 
Using the number of standard deviations is therefore a more ro- 
bust approach and is consistent with our previous work (M cEwen 
et al. |2005a| ). Significance levels corresponding to the detections 
considered in Fig. ^ are calculated and displayed in Table ^ For 
clarity, we show only those values from the co-added WMAP map, 
although the ILC map exhibits similar results. These results con- 
firm our inferences from direct observation of the statistics relative 
to the confidence levels and histograms shown in Fig. ^ and Fig. ^ 
respectively: the original kurtosis detections of non-Gaussianity are 
eliminated, while the original skewness detections remain. We also 
determine the significance of the new candidate detections of non- 
Gaussianity observed in the kurtosis of the Mexican hat wavelet 
coefficients in the Bianchi corrected data. Of the 1000 simual- 
tions, 115 contain a statistic that exhibits a greater deviation that 
the symmetric Mexican hat wavelet kurtosis on scale ag of the 
Bianchi corrected data, hence this detection may only be made at 
the 88.5% significance level. 448 of the simulations contain a statis- 
tic that exhibits a greater deviation that the elliptical Mexican hat 
wavelet kurtosis on scale au of the Bianchi corrected data, hence 
this candidate detection may only be made at the 55.2% signifi- 
cance level. Thus we conclude that no highly significant detection 
of non-Gaussianity can be made on any scale or orientation from 



the analysis of the kurtosis of spherical wavelet coefficients in the 
Bianchi corrected data, however the detections made previously in 
the skewness of spherical wavelet coefficients remain essentially 
unaltered. 

Finally, we perform tests to probe the significance of de- 
viations from Gaussianity in the aggregate set of test statistics. 
These tests inherently take all statistics on all scales and orien- 
tations into account. The results of these tests are summarised in 
Fig. ^ The overall significance of the detection of non-Gaussianity 
is reduced for the Mexican hat wavelets, although this reduction is 
not as marked as that illustrated in Table ^ since both skewness and 
kurtosis statistics are considered when computing the;^^^, and it is 
only the kurtosis detection that is eliminated. For example, when an 
equivalent test is performed using only the kurtosis statistics the 
significance of the detection made with the symmetric Mexican hat 
wavelet drops from 99.9% to 95% (note that this is still consider- 
ably higher than the level found with the previous test, illustrating 
just how conservative the previous method is). The significance of 
the detection made with the real Morlet wavelet is not affected by 
correcting for the Bianchi template. This is expected since the de- 
tection was made only in the skewness of the real Morlet wavelet 
coefficients and not the kurtosis. 

We quote the overall significance of our detections of non- 
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Table 2. Deviation and significance levels of spherical wavelet coefficient 
statistics calculated from the WMAP and Bianchi corrected WMAP maps 
(similar results are obtained using the ILC map). Standard deviations and 
significant levels are calculated from 1000 Monte Carlo simulations. The 
table variables are defined as follows: the number of standard deviations 
the observation deviates from the mean is given by No-', the number of sim- 
ulated Gaussian maps that exhibit an equivalent or greater deviation in any 
test statistics calculated using the given wavelet is given by N^q^\ the cor- 
responding significance level of the non-Gaussianity detection is given by 
6. Only those scales and orientations corresponding to the most significant 
deviations from Gaussianity are listed for each wavelet. 



Skewness 
(«2 = 1000 



Kurtosis 
(«6 = 3000 





WMAP 


WMAPb.viIh 


WMAP 


WMAPB.viih 


A^dev 

6 


-3.38 
28 maps 
97.2% 


-3.53 
21 maps 
97.9% 


3.12 
47 maps 
95.3% 


1.70 
605 maps 
39.5% 






(a) Mexican hat e 


= 0.00 






Skewness 
(fl3 = 150'; 7 = 72°) 
WMAP WMAPb.viIh 


Kurtosis 
(flio = 500';7= 108°) 
WMAP WMAPeviih 


Na 

6 


-4.10 

39 maps 
96.1% 


-4.25 
29 maps 
97.1% 


3.01 
199 maps 
80.1% 


1.88 
887 maps 
11.3% 






(b) Mexican hat e 


= 0.95 






Skewness 
(an = 550'; 7 = 72°) 
WMAP WMAPB.viih 


Kurtosis 
(an = 550'; 7 = 72°) 
WMAP WMAPB.viih 


Na 
A^dev 

5 


-5.61 
17 maps 
98.3% 


-5.66 
16 maps 
98.4% 


2.66 
642 maps 
35.8% 


2.67 
628 maps 
37.2% 



(c) Real Morlet A: = (10, 0)^ 



Gaussianity at the level calculated by the first approach, since this 
is the more conservative of the two tests. 




WMAP 
((5 = 99. £ 



((5 = 99.7%) 



ILC 

((5 = 99.9%) 



(a) Mexican hat e = 0.00 




X / DOF 

(b) Mexican hat e = 0.95 



Li 



VI MAP 
(6 = 99.8%) 



= 99.9%) 
VIlMAP- BIANCHI 



-BIANCH 
3.9%) 



3.3 Localised deviations from Gaussianity 

The spatial localisation inherent in the wavelet analysis allows one 
to localise most likely sources of non-Gaussianity on the sky. We 
examine spherical wavelet coefficients maps thresholded so that 
those coefficients below 3cr (in absolute value) are set to zero. The 
remaining coefficients show likely regions that contribute to devia- 
tions from Gaussianity in the map. 

The localised regions of the skewness-flagged maps for each 
wavelet are almost identical for all of the original and Bianchi cor- 
rected co-added WMAP and ILC maps. This is expected as it has 
been shown that the Bianchi correction does not remove the skew- 
ness detection. All of these thresholded coefficient maps are almost 



(c) Real Morlet A: = (10, 0)^ 



Figure 5. Histograms of normahsed x test statistics obtained from 1000 
Monte Carlo simulations. Normalised x'^ values corresponding to the fol- 
lowing maps are also plotted: WMAP combined map (solid, blue, square); 
ILC map (solid, green, circle); Bianchi corrected WMAP combined map 
(dashed, blue, triangle); Bianchi corrected ILC map (dashed, green, dia- 
mond). The significance level of each detection made using x'^ values is 
also quoted (S). 
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identical to those shown in Fig. 9 (a,c,d) of McEwen et al. (2005a), 
hence they are not shown here. 

The locaUsed regions detected in the kurtosis-flasged maps 
for the Mexican hat wavelets are shown in Fig. (the real 
Morlet wavelet did not flag a significant kurtosis detection of 
non-Gaussianity). The thresholded coefficient maps for the orig- 
inal and Bianchi corrected data are reasonably similar, however 
the size and magnitude of the cold spot at Galactic coordinates 
(l,b) = (209° , -57°) is s i gnific antly reduced in the Bianchi cor- 
rected maps. |Cruz et aL (2005) claim that it is this cold spot that 
is responsible for the kurtosis detection of non-Gaussianity. This 
may explain why the kurtosis detection of non-Gaussianity is elim- 
inated in the Bianchi corrected maps. 

Although the new cadidate detections of non-Gaussianity ob- 
served in the kurtosis of the Mexican hat wavelet coefficients of the 
Bianchi corrected data were shown in section 3^ not to be particu- 
larly significant, we nevertheless construct the localised maps cor- 
responding to these candidate detections. The regions localised in 
these maps show no additional structure than that shown in Fig. |^ 
The only significant difference between the localised regions is that 
the cold spot at (/, b) = (209°, -57°) is absent. 




(a) Gaussian plus Bianchi simulated 
map Mexican hat wavelet coefficients 
(6 = 0.00; ae = 3000 




3.4 Gaussian plus Bianchi simulated CMB map 

It addition to testing the WMAP data and Bianchi corrected ver- 
sions of the data, we also consider a simulated map comprised of 
Gaussian CMB ffuctuations plus an embedded Bianchi component. 
We use the same strategy to simulate the Gaussian component of 
the map that is used to create the Gaussian CMB realisations used 
in the Monte Carlo analysis and add to it a scaled version of the 
Bianchi template that was fitted to the WMAP data by J05. The mo- 
tivation is to see whether any localised regions in the map that con- 
tribute most strongly to non-Gaussianity coincide with any struc- 
ture of the Bianchi template. 

Non-Gaussianity is detected at approximately the 3cr level in 
the kurtosis of the symmetric and elliptical Mexican hat wavelet 
coefficients once the amplitude of the added Bianchi template is 
increased to approximately (g)^ ~ 15 x 10~^^ (approximately four 
times the level of the Bianchi template fitted by J05), corresponding 
to a vorticity of ~ 39 X 10~^^. No detections are made in any 
skewness statistics or with the real Morlet wavelet. The localised 
regions of the wavelet coefficient maps for which non-Gaussianity 
detections are made are shown in Fig. ^. The Mexican hat wavelets 
extract the intense regions near the centre of the Bianchi spiral, with 
the symmetric Mexican hat wavelet extracting the symmetric struc- 
ture and the elliptical Mexican hat wavelet extracting the oriented 
structure. This experiment highlights the sensitivity to any Bianchi 
component of the Mexican hat kurtosis statistics, and also the in- 
sensitivity of the Mexican hat skewness statistics and real Morlet 
wavelet statistics. The high amplitude of the Bianchi component re- 
quired to make a detection of non-Gaussianity suggests that some 
other source of non-Gaussianity may be present in the WMAP data, 
such as the cold spot at (l,b) = (209°, -57°), and that the Bianchi 
correction may act just to reduce this component. 

3.5 Foregrounds and systematics 

From the proceeding analysis it would appear that a Bianchi com- 
ponent is not responsible for the non-Gaussianity observed in the 
skewness of the spherical real Morlet wavelet coefficients. The 
question therefore remains: what is the source of this non-Gaussian 



(b) Gaussian plus Bianchi simulated 
map Mexican hat wavelet coefficients 
(e = 0.95;flio = 500'; 7= 108°) 

Figure 7. Thresholded Mexican hat wavelet coefficients of the Guassian 
plus Bianchi simulated map. The coefficient maps shown are flagged by a 
kurtosis detection of non-Gaussianity. Notice how the Mexican hat wavelets 
extract the intense regions near the centre of the Bianchi spiral, with the 
symmetric Mexican hat wavelet (e = 0.00) extracting the symmetric struc- 
ture and the elliptical Mexican hat wavelet (e = 0.95) extracting the oriented 
structure. 



signal? We perform here a preliminary analysis to test whether un- 
removed foregrounds or WMAP systematics are responsible for the 
non-Gaussianity. 

The coadded map analysed previously is constructed from a 
noise weighted sum of two Q-band maps observed at 40.7GHz, 
two V-band maps observed at 60.8GHz and four W-band maps ob- 
served at 93.5GHz. To test for foregrounds or systematics we ex- 
amine the skewness observed in the separate WMAP bands, and 
also in difference maps constructed from the individual bands. In 
Fig. js] (a) the skewness of real Morlet wavelet coefficients is shown 
for the individual band maps Q = Q1+Q1,V = V1+V2 and 
W = Wl + W2 -h W3 + W4, and in Fig. js] (b) the skewness is shown 
for the diff^erence maps VI - V2, Ql - Q2, Wl - W4 and W2 - W3 
(the W-band diff'erence maps have been chosen in this order to en- 
sure that the beams of the maps compared are similar). Note that the 
conffdence regions shown in Fig. (b) correspond to the WMAP 
coadded map and not the diff'erence maps. It is computationally 
expensive to compute simulations and signiffcance regions for the 
diff'erence maps, thus one should only compare the skewness sig- 
nal with that observed previously. One would expect any detection 
of non-Gaussianity due to unremoved foregrounds to be frequency 
dependent. The skewness signal we detect on scale an is identi- 
cal in all of the individual WMAP bands, hence it seems unlikely 
that foregrounds are responsible for the signal. Moreover, since the 
skewness signal is present in all of the individual bands it would 
appear that the signal is not due to systematics present in a single 
WMAP channel. The signal is also absent in the difference maps 



8 /. D. McEwen et al. 




(a) Co-added WMAP map Mexican hat wavelet coefficients (e = 
0.00; ae = 3000 




(c) Co-added WMAP map Mexican hat wavelet coefficients (e = 
0.95; flio = 500^ 7 = 108°) 




(b) Bianchi corrected co-added WMAP map Mexican hat wavelet co- 
efficients (e = 0.00; ae = 3000 




(d) Bianchi corrected co-added WMAP map Mexican hat wavelet co- 
efficients (6 = 0.95; aio = 500'; y = 108°) 



Figure 6. Thresholded spherical wavelet coefficients for the original and Bianchi corrected co-added WMAP map. The inset figure in each panel shows a 
zoomed section (of equivalent size) around the cold spot at (/, b) = (209°, -57°). The size and magnitude of this cold spot is reduced in the Bianchi corrected 
data. Only those coefficient maps corresponding to the most significant kurtosis detections fo r the Mexican hat wavelets are shown. Other coefficient maps 
show no additional information to that presented in our previous work ( [vicEwen et al. I I2OO54 (see text). The corresponding wavelet coefficient maps for the 
ILC map are not shown since they are almost identical to the coefficients of the co-added WMAP maps shown above. 



that are dominated by systematics and should be essentailly absent 
of CMB and foregrounds. 

From this preliminary analysis we may conclude that it is un- 
likely that foregound contamination or WMAP systematics are re- 
sponsible for the highly significant non-Gaussianity detected with 
the spherical real Morlet wavelet. This analysis has also highlighted 
a possible systematic in the Q-band on scale a^. A more detailed 
analysis of this possible systematic and a deeper analysis of the 
cause of the non-Gaussianity detected with the real Morlet wavelet 
is left for a separate piece of work. 



4 CONCLUSIONS 

We have investigated the eflfect of correcting the WMAP data for 
a Bianchi type Vllh template on our previous detections of non- 
Gau ssianity made with directional spherical wavelets (M cEwen 
et al. [2005a| ). The best-fit Bianchi template was simulated with the 
parameters determined by JOS using the latest shear and vorticity 
estimates (J05; private communication). We subsequently used this 
best-fit Bianchi template to 'correct' the WMAP data, and then re- 
peated our wavelet analysis to probe for deviations from Gaussian- 
ity in the corrected data. 

The deviations from Gaussianity observed in the kurtosis of 
spherical wavelet coefificients disappears after correcting for the 
Bianchi component, whereas the deviations from Gaussianity ob- 
served in skewness statistics are not afifected. The highly signifi- 
cant detection of non-Gaussianity previously made in the skewness 
of real Morlet wavelet coefficients remains unchanged at 98% us- 
ing the extremely c onse rvative method to compute the significance 
outlined in section 3.2, The tests also performed indicate that 



the Bianchi corrected data still deviates from Gaussianity when all 
test statistics are considered in aggregate. Since only the skewness- 
flagged detections of non-Gaussianity made with the Mexican hat 
wavelet remain, but the kurtosis ones are removed, the overall sig- 
nificance of Mexican hat wavelet tests are reduced. There was 
no original detection of kurtosis in the real Morlet wavelet coeflfi- 
cients, thus the significance of the remains unchanged for this 
wavelet. Finally, note that one would expect the skewness statistics 
to remain unaffected by a Bianchi component (or equivalently the 
removal of such a component) since the distribution of the pixel 
values of a Bianchi component is itself not skewed, whereas a sim- 
ilar statement cannot be made for the kurtosis. 

Regions that contribute most strongly to the non-Gaussianity 
detections have been localised. The skewness-flagged regions of 
the Bianchi corrected data d o not diflfer sig n ificantl y from those 



regions previously found in 



McEwen et al 



( |2005a| ). 



One would 



expect this result: if these regions are indeed the source of non- 
Gaussianity, and the non-Gaussianity is not removed, then when 
most likely contributions to non-Gaussianity are again localised the 
regions should remain. The kurtosis-flagged regions localised with 
the Mexican hat wavelets are not markedly altered by correcting for 
the Bianchi template, however the size and magnitude of the cold 
spot at Galactic coordina tes (/, b) = (209°, -57°) is significantly re- 
duced. |Cruz et al" (2005) claim that it is solely this cold spot that 
is responsible for the kurtosis detections of non-Gaussianity made 
with Mexican hat wavelets, thus the reduction of this cold spot 
when correcting for the Bianchi template may explain the elimi- 
nation of kurtosis in the Bianchi corrected maps. 

After correcting the WMAP data for the best-fit Bianchi Vllh 
template, the data still exhibits significant deviations from Gaus- 
sianity, as highlighted by the skewness of spherical wavelet coef- 
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(a) Individual WMAP band maps: WMAP co- 
added (solid, blue, square); Q (solid, green, circle), 
V (dashed, blue, triangle); W (dashed, green dia- 
mond). 




(b) WMAP band difference maps: Ql - Q2 (solid, 
blue, square); VI - V2 (solid, green, circle); Wl - 
W4 (dashed, blue, triangle); W2 - W3 (dashed, 
green, diamond). 

Figure 8. Skewness for individual and difference WMAP band maps - real 
Morlet k = (10,0)^. Note that the strong non-Gaussianity detection made 
on scale is present in all of the individual band maps but is absent from 
all of the difference maps that should contain predominantly systematics. 
The confidence regions shown in these plots are for the WMAP coadded 
map (see comment text). 



ficients. A preliminary analysis of foreground contamination and 
WMAP systematics indicates that these factors are also not respon- 
sible for the non-Gaussianity. A deeper investigation into the source 
of the non-Gaussianity detected is required to ascertain whether the 
signal is of cosmological origin, in which case it would provide ev- 
idence for non-standard cosmological models. Bianchi models that 
exhibit a small universal shear and rotation are an important, alter- 
native cosmology that warrant investigation and, as we have seen, 
can account for some detections of non-Gaussian signals. How- 
ever, the current analysis is only phenomenological since revisions 
are required t o update the B i anchi -induced temperature fluctuations 
calculated by iBarrow et al.| ( |l985| ) for a more modern setting. Nev- 
ertheless, such an analysis constitutes the necessary first steps to- 
wards examining and raising the awareness of anisotropic cosmo- 
logical models. 
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APPENDIX A: BIANCHI TYPE Vllh REPRESENTATIONS 

We do not attempt to describe or solve the geodesic equations re- 
sulting from t he Bi anchi Vllh model, but rather refer to reader to 



Barrow et alj ( |1985| ) for details. Enough detail on the induced tem- 
perature fluctuations is presented here so that the interested reader 
may re-produce our simulated maps. 

Analytic forms are given to compute the Bianchi-induced tem- 
perature fluctuations directly in either real or harmo nic space. We 



correct some errors (most likely typographical) in Barrow et al 



(1985) in the analytic form of the harmonic space representation. 



Different representations are preferred for different applications, 
hence the ability to compute directly in the space required is ob- 
viously beneficial. The efficiency of computing in either space is 
not markedly different, however computing the harmonic coeffi- 
cients directly does have some advantages. For instance, the 6 reso- 
lution requi red to compute numerically some integrals (specifically, 
(A18) and (|a19|)) is independent of the resolution of the real space 



map, and typically the integrals may be computed accurately for 
a relatively low 6 resolution. Moreover, the rotation may be per- 
formed considerably more efficiently and accurately^ in harmonic 
space. Finally, the Bianchi-induced fluctuations have a low band- 
limit, thus harmonic coefficients need only be computed for a rela- 
tively low Cax- 



Al Real space representation 

The Bianchi Vllh -induced temperature fluctuations are parame- 
terised by the total energy density ^^o, the redshift of recombina- 



^ A finite point set on the sphere that is invariant under rotations does not 
exist, thus arbitrary rotations cannot be performed exactly in real space on 
a pixelised sphere. 



tion ze and the Bianchi shear (f)^, handedness k and x parame- 
ters, where H is the Hubble parameter. The handedness parameter 
takes the values k - +\ and k - -\ fox maps with right and left 
handed spirals respectively. Physically, x is related to the character- 
istic wavelength over which the principle axes of shear and rotation 
change orientation; hence x defines the 'spiralness' of the resul- 
tant temperature fluctuations. The h parameter of the Bianchi Vllh 
model is related to x and Qq by 



The vorticity of a Bianchi Vllh model is related to the other param- 
eters by 



a- 



V2(l + hyi^il + 9h) 



1/2 



(A2) 



6x2Qo 

The temperature fluctuations induced in the Bianchi type Vllh 
model are given by 



-{Oo, 0o) = (^)^ {[A(^o) + BiOo)] sin(0o) 
+ /c[m)-A(^o)] cos(0o)}, 



(A3) 



where and 0o specify the direction of the incoming photon, 
which is related to the observing angles by 6o\y - n - Oq and 
0ob = + 00- M^o) and B(6o) are defined by 

A(6o) = Ci sin(^o) 

- CiC2(^o) {cos [i//(te, Oo)] - 3/1^/2 sin [i//(te, ^o)]} 
-0 s(t, 6o) [l - s\t, 6o)] sin [i//(t, 6o)] dr 



and 



+ 0. 



?>h^'^Ci sin(^o) 



TE [1 + ^2(t, 6*0)] sinh^(/zi/2T/2) 



(A4) 



- CiC2(6'o) {sin [j/r(T£, ^o)] + 3/^^/2 cos [^/^(t^, 6*0)]} 
-0 5(t, ^0) [1 - s\t, Oo)] cos [jA(t, ^0)] dr 



JTf 



TE [1 + ^2(t, 6*0)] sinh^(/zi/2T/2) 



,(A5) 



where we integrate over the photon path in conformal time r from 
photon emission (last scattering surface) 



te = 2h-^l^ sinh 



1 +Ze 



1/2 



xl/2 



to observation 

To = 2r^/2 sinh"^ [{QT^^ - If ^ 
The terms Ci, C2(^o)» C'b are defined by 
Ci = (3^^o-^)"' , 



C2(^0) ^ 



2s{TE,eo){\+ZE) 



l + s\TE.eQ) ' 

and 

C3 = 4/^^/2(1 - ^of'^or^^ . 

The functions s{t, 6q) and j/^(t, ^0) are defined by 



tan(^o/2) 



^(T,^o) 
and 

V^(T,^0) = (T-To)-r 

, _2/7l/2(T-T0) 



1/2 



In [sin2(^o/2) 
s'(^o/2)] . 



(A6) 

(A7) 

(A8) 
(A9) 

(AlO) 
(All) 

(A12) 
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At a sacrifice to notational simplicity, but, we hope, to add clarity, 
we have made the dependence on all variables that are not intrinsic 
Bianchi parameters explicit in all functions. 

The induced swirl pattern described by these equations is cen- 
tred on the south pole. Obviously, the centre of the swirl may by 
located at any position in our coordinate system, thus a final ro- 
tation is required to fully describe any possible realisation of the 
Bianchi-induced temperature fluctuations. This introduces three ad- 
ditional parameters, the three Euler angles that describe the rotation 
{a,p, y) (for which we adopt the active zjz-convention, where func- 
tions rather than axes are rotated). 

A2 Harmonic space representation 

We derive the analytic form of the spherical harmonic coefficients 
of the Bianchi-induced temperature fluctuations in this section. Due 
to the low frequency structure of the Bianchi temperature fluctua- 
tions, both globally and azimuthally, the band-limit of harmonic 
coefficients is low in both i and m. Thus, computing simulated 
temperature fluctuations may be performed efficiently and, as we 
explain, more accurately in harmonic space. 

We decompose the Bianchi-induced temperature fluctuations 
into a sum of spherical harmonics 

— (^, 0) = ^ ^ a,^YUO. 0) , (A13) 
^ e=\ m=-e 

where the harmonic coefficients are given by the usual projection 
on to each spherical harmonic basis function 



(A14) 



and dn = sin 6 dO d0 is the usual rotation invariant measure on the 
two-sphere S^.ln practice the outer summation of ( A13 ) is trun- 
cated to ^max terms. We adopt the Condon-Shortley phase conven- 
tion, where the normalised spherical harmonic are defined by 



Jo 
and 



P-"\X) 



*cos0d0 = n5m,±\ 



(^ + m)! 



(A21) 



(A22) 



Notice that the minimal azimuthal structure of the induced tempera- 
ture fluctuations ensures that the harmonic coefficients are non-zero 
only for m = +1. For completeness, we also state the analytic form 
for the power spectrum of the Bianchi-induced fluctuations: 

The swirl pattern of the induced temperature fluctuations is 
again centred on the south pole in the expression given by ( |\17| ). 
We may perform the required rotation directly in harmonic space, 
noting that the spherical harmonics of a rotat ed function a re related 
to the harmonics of the original function by ( ^isbc||l996| ) 



(A24) 



where , are the Wiener D-matrices. The harmonic coefficients 

m,m c5 

of the rotated Bianchi-induced temperature fluctuation map simply 
reduce to 



(A25) 



since the harmonic coefficients of the original map are non-zero 
only for m - ±1. 



r,.(^,0) = (-ir 



2^+1 {l-m)\ 
47T {£ + m)\ 



(A15) 



where P'^(x) are the associated Legendre functions. If we adopt the 
photon, rather than observing angles for the temperature fluctua- 
tions we simply introduce a factor of (-1)^+'": 



- I — 

J52 ^0 



(A16) 



Substituting t he ex pression for the temperature fluctuations given 
by into (|a14|) yields 



a„n = (^)^ 7:[mK(lf - /f ) + + /f )] V±i 
where Sij is the Kronecker delta. 



r - 



and 



\4;r^(^+l) Jo 



/ 2^ + 1 r 
\ WTT) Jo 



A(6')Pj (cos 6*) sin6'd6' 



B(e)Pl(cose) sin6'd6'. 



(A17) 



(A18) 



(A19) 



In the derivation of (A17) we have made use of the following rela- 
tions: 



f 

Jo 



-im(j) 



sin0d0 - -im7T6m,±i 



(A20) 



